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ly-j | Abstract. Two XMM-Newton observations of the black-hole binary GRS 1915+105 were triggered in 2004 (April 17 and 21), 

\Q ■ during a long "plateau" state of the source. We analyzed the data collected with EPIC-pn in timing and burst modes, respectively. 

Reflection Grating Spectrometers were used only on April 21st. The source 2-10 keV flux is ~ 0.6 (unabsorbed: 0.9 -f 1.1) 
xl0~ 8 in cgs units. While the light curves show only small amplitude variations (a few percent) at timescales longer than a few 
seconds, a QPO is seen at about 0.6 Hz - as expected in,^ variability modes of GRS 1915+105, when the phenomenological 
correlation with the source flux is taken into account - possibly with a harmonic signal at 1.2 Hz. The pn spectrum is well fitted 
without invoking thermal disk emission, on the basis of four main components: a primary one (either a simple power law or 
O j 1 thermal Comptonization models), absorbed by cold matter with abundances different than those of standard ISM; reprocessing 

from an ionized disk; emission and absorption lines; and a soft X-ray excess at ~ 1 keV. However, the last is not confirmed by 
the RGS spectra, whose difference from the EPIC-pn ones lacks a fully satisfactory explanation. If real, the soft X-ray excess 
may be due to reflection from an optically thin, photoionized disk wind; in this case it may lead to a way to disentangle intrinsic 
. from interstellar absorption. 

Key words. Black hole physics - Accretion, accretion disks - X-rays: binaries - X-rays: individuals: GRS 1915+105 
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. 5r , 1. Introduction monly adopted value of the source's inclination (70°; but see 

Maccarone, 2002, for a discussion), Greiner et al. (2001b) used 

The well-known black-hole (BH) binary GRS 1915+105 is thesa me observations to determine the mass of the central com- 

somettmes classified as a superluminal mwroquasar because pact objecti which has been constrained t0 Mc = 14 + 4M Q , 

of the impressive radio jets whose matter seems to move with i e weU aboye ±e maximum neutron star mass Umit Thus 

superluminal velocity due to a relativistic effect (Mirabel & GRS 1915+105 is believed t0 host a BH with a grav i tationa i 

Rodriguez, 1994; Rodriguez & Mirabel, 1999; Fender et al., radiuf , , g = cm ^ 2 1 km, so it is classified as a BH binary. 

1999). For a recent review of this source, including a general . , c ,,,,,, 

,. . „. . „ , „ „ . ,r,nrm\ Reig et al. (2003) note that, although the phenomenology of 

discussion of its properties, see Fender & Belloni (2004). ^^o,^,,^ „ i , . , 

GRS 1915+105 has been usually described in terms of three 

Due to very large obscuration, the spectral type of . , . . ■ . D , „ , „ , • , , 

3 & 3V spectral states named A, B, and C (the latter being low/hard 

GRS 1915+105's companion (a K-M III star) was discovered , . . , . . ,x . t . , .• , 

F v ' and jet-dominated), in fact several source properties always 

only recently, via infrared observations (Greiner et al., 2001). , . • , . . . / .• „ ■ 

3 ■" ' correspond to the canonical intermediate (sometimes called 

This yielded evidence that the system, initially discovered in , . , x , . c , .. nTT - TT r> n • 

3 very high) state of galactic BH sources (see Homan & Belloni, 

the y rays (Castro-Tirado et al. 1992), belongs to the class of ^ and Belloni ^ ^ ^ for a definition) This is con . 

low-mass X-ray binaries (LMXBs). On the basis of the com- c ■ , ^ . , / OAA /i\ u i t .l tU > 
3 ' firmed by Done et al. (2004), who also note that the source s 

unique limit-cycle variability appears when the source radiates 

at super-Eddington luminosities. Many of the GRS 1915+105 
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by far the longest of any LMXB (Greiner et al., 2001) indi- 
cate that the Roche lobe overflow occurs in a very wide binary. 
The cause of all the unique long, as well as short-term, vari- 
ability of GRS 1915+105 is then linked to the evolution of the 
huge disk structure, which can contain enough material to man- 
tain nearly-Eddington accretion rates over timescales of several 
years. 

A strong and broad iron Ka line was present in BeppoSAX 
(Martocchia et al. 2002, 2004) and ASCA GIS (Miller et al. 
2004) spectra of GRS 1915+105. Rossi-XTE PC A data are 
generally well-fitted including a broad iron line component 
(e.g. Done et al. 2004), but the energy resolution of this instru- 
ment does not allow discriminating among different line mod- 
els. The Fe Ka emission was also compatible with the spec- 
tral fits of Chandra data, taken during a low/hard state of the 
source, as reported by Lee et al. (2002). 

An XMM-Newton ToO observation of GRS 1915+105, 
which would allow us to study the iron line profile with en- 
hanced spectral resolution and sensitivity, was proposed in the 
second Announcement of Opportunity (A02). The observation 
was intended to be triggered by the occurrence of a "plateau" 
state of the source (Foster et al. 1996), similar to the one ob- 
served during the BeppoSAX 1998 observation; this was also 
necessary in order to have the source in a less dramatic variabil- 
ity state and at a lower flux level to minimize technical prob- 
lems due to instrumental pile-up and telemetry. For the latter 
reason we also planned to observe in timing mode. 

The observation was triggered in April 2004, divided into 
two parts with ~ 4 days separation one from the other. In this 
paper, while presenting the results, we offer an XMM-Newton 
view of GRS 1915+105 for the first time. First of all we de- 
scribe (in Sect. 2) the observations and data extraction meth- 
ods, and details of the state of the source based on its variabil- 
ity and general considerations. In Sect. 3 the best-fit spectral 
models and related physical parameters are given. Finally, in 
the conclusions, we summarize these XMM-Newton measure- 
ments and put our results in the context of the known source 
properties and open issues. 

2. Observations and data reduction 

To trigger the ToO observation, we set a condition on variabil- 
ity (rms single-dwell points of Rossi-XTE ASM to stay < 9% 
within one day) and a condition on flux (Rossi-XTE ASM daily 
rate below a threshold of 30 cts/s). 2 Both triggering conditions 
were met during the long "plateau", a lower-flux state also 
characterized by a more intense radio emission (Foster et al. 
1996), which lasted several weeks during the Spring of 2004 
(Fig. 1). The observation could be therefore started, in EPIC 
timing mode, on April 17, 2004 (hereafter: OBS1). 

However, the source counts were enough to saturate the 
telemetry, even after switching off the MOS and RGS. Thus 
a second part of the observation was performed on April 21, 
2004 (hereafter: OBS2), this time in (EPIC-pn) burst mode. 

2 ASM is the All Sky Monitor onboard Rossi-XTE. We used 
the public ASM "Weather Map" mask, accessible at: http:// 
heasarc.gsfc.nasa.gov/xte_weather/. 
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Fig. 1. Radio and X-ray light curves of GRS 1915+105 before and dur- 
ing the Spring 2004 plateau (starting at ~ 2453070 JD), from: http:// 
www.mrao.cam.ac.uk/~guy/. Radio data at 15 GHz were observed 
with the Ryle Telescope, Cambridge (see Pooley & Fender, 1997, for 
a description of the techniques used); X-ray data (ASM flux and hard- 
ness ratio) are taken from the Rossi-XTE quick-look web pages: http:// 
xte.mit.edu/ASM_lc.html. The approximate dates of OBS1 and OBS2 
are shaded. Note the anti-correlation between X and radio. 



Given the typical plateau behavior, with strong jet activity 
and high radio level, we conclude that the source was caught 
in the conventional "hard" (C) spectral state and x variability 
class, as defined by Belloni et al. (1997a,b, 2000), in both ob- 
servations. This is similar to the condition in which the source 
was seen by simultaneous INTEGRAL (rev.57) and Rossi-XTE 
observations performed in the Spring of 2003 (Hannikainen 
et al. 2003, 2004, 2005; see Fuchs et al. 2003 for a multi- 
wavelength study of one of these observations). 3 

Count rate thresholds for pileup in timing modes are re- 
ported in the XMM User's Handbook. 4 The MOS count rate 
registered in both observations (~ 130 and 160 cts/s) is well 
above the quoted limit (100); therefore EPIC-MOS data are 
unusable, 5 and we restricted the analysis to EPIC-pn and RGS 
data (the latter are available only for OBS2). 

Since in this state the source is stable, we can use the whole 
integration time to obtain very high statistics and signal/noise. 
Through adopting the customary method of extracting high- 
energy (E > 10 keV) background light curves, we checked for 
intervals of low contamination in both observations. While no 
background flaring is apparent during OBS1, we did initially 
select a 9 ks low-background interval (198951000+ 198960000 
s) in OBS2. We performed preliminary OBS2 data analysis 
with this selected dataset, then afterwards verified that the re- 
sults were not modified when considering the whole observ- 

3 This radio-loud state had been described already by Muno et al. 
2001, and by Trudolyubov 2001 (as "type II"). 

4 See Table 3 in: Ehle et al. (2004) - http:// 
xmm.vilspa.esa.es/external/xmm_user_support/documentation/ 
uhb/node28.html#2828 

5 Excising the core of the PSF is enough to cope with this problem 
in the case of imaging modes, but it is not applicable for timing modes. 
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ing interval. Therefore, in the following, we report result 
the data analysis that was performed using all observing ti 
e.g. with the highest possible statistics. This also allowec 
to refine the spectro-variability analysis (e.g. the analysi: 8 
root mean squared variability vs. energy: Ponti et al. 2( 
as much as possible, which we performed to possibly dif 
tangle the various contributions to the observed spectra ; 
thus, to get an independent validation of the adopted er 
sion/absorption model. The high number of counts allow?. § 
high-resolution analysis of rms vs. E over the whole EF " 
pn energy band (about 35 energy bins for OBS2, and sev 
dozen energy bins for OBS 1 ; the minimum timescales of va 
tion we could consider were, respectively, ~ 1000 and ~ 10( ° 
We used background-subtracted light curves in order to ai 
contamination at the lowest and highest energies (E < 2 
E > 7 keV), and found a mean rms value on the order of 
(0.02 4- 0.07) in both observations, compatible with a null 
when Poisson noise is removed, with no visible features in 
"rms spectra". We further verified this by studying F var vs. 

The data were reduced with the standard SAS v6.1 and 
FTOOLS software packages. For EPIC-pn, single and dou- 
ble events were used in the timing modes (PATTERN <= 4); 
the fit was restricted to energies > 0.5 (> 0.4) keV in timing 
(burst) mode in order to avoid the increased noise. As recom- 
mended by the XMM-Newton Helpdesk, no selection on the 
RAWY coordinate was applied since this is related to a fine- 
time and the selection would therefore exclude certain time pe- 
riods. However, RAWY < 160 was imposed as a condition in 
burst mode to avoid direct illumination by the source. Response 
and auxiliary matrices were created with the SAS tools rmfgen, 
arfgen, rgsrmfgen. 

For the fit procedure we assumed a source distance of 12.5 
kpc and inclination ; ~ 70° (Rodriguez & Mirabel, 1999; 
Fender et al., 1999). 

2.1. First observation 

The April 17 observation (OBS1: XMM-Newton rev.798, 
start time: 14:18:56 [JD24531 13.096], stop time: 20:03:37 
[JD24531 13.336]) was performed in (both pn and MOS1) tim- 
ing mode, as requested in the A02 proposal mainly with the 
aim of avoiding pile-up problems. Moreover, the MOS2 cam- 
era and both RGS had to be switched off, to avoid overloading 
the telemetry. Nevertheless, the source flux was still so high 
that telemetry drops were frequent even in pn data, resulting in 
short gaps in the light curves, lasting a few to several seconds, 
for a total of about 75% of the observing time. They cause the 
search for (quasi-) periodic oscillations to be senseless on ~ 1 
Hz timescales. Source light curves do not show any apprecia- 



GRS 1915+105 (April 21. 2004) 



Bin time: 40.96 ! 



6 The nominal time resolution is 29 and 7 microseconds for EPIC- 
pn timing and burst modes, respectively. However, all results derived 
from the light curves must be taken with caution: in the first obser- 
vation, even after removing all telemetry drops, the collected photons 
"saturated" the instrument's bandwidth; therefore, energy-dependent 
oscillations in the count rate may have been "cut off". On the other 
hand, photons collected in the second observation, i.e. in burst mode, 
are selected on the basis of an ad hoc duty cycle in the reading process 
(events are not registered for 97% of the time). 




Fig. 2. The source light curve in the 0.4 - 10.0 keV band during OBS2 
(EPIC-pn, burst mode, not background-subtracted). The whole 20946 
s observation has a mean 618 (i.e. 18.54/0.03) cts/s; the background 
has just about 18 (i.e. 0.5320/0.03) cts/s. 



ble variations at timescales ~ a few to ~ 100 s; however, as ex- 
plained above (see Footnote 6), due to the full scientific buffer 
and telemetry drops, this result should be taken with caution. 

Even though affected by telemetry drops for up to about 
75% of the observing time, EPIC-pn data are well below the 
pile-up count rate threshold 7 and can therefore be safely used 
for spectral analysis. 

The data were selected from a four-column strip centered 
on the source throughout the timing-mode (one-dimensional) 
"image"; a similar strip was selected far away from the source 
as a background sample. 

2.2. Second observation 

While the source essentially persisted in the same state, a sec- 
ond observation was started on April 21st (OBS2, rev. 800, 
start time: 13:49:23 [JD24531 17.076], stop time: 20:56:55 
[JD24531 17.373]). To cope with the telemetry problem, this 
time we used EPIC-pn in burst mode, MOS in timing mode, 
and RGS in high count rate (HCR) mode. The other cameras 
were off. 

In this observation, a background flare is apparent at the 
very beginning both in EPIC and RGS data; however, due to 
the very high S/N and shortness of the flare, it does not re- 
ally affect the spectra. Source data were selected from a six- 
column strip (i.e. two columns more than in timing mode, in 
order to increase the statistics) throughout the burst-mode (one- 
dimensional) "image"; a similar strip was selected, far away 
from the source, as a background sample. 



7 For EPIC-pn timing mode the pile-up threshold for a point source 
is 1500 cts/s, the telemetry threshold is 300 cts/s, while we registered 
~ 650 cts/s. 
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Fig. 3. Power spectrum in the 0.4 - 13.0 keV band of the OBS2 (April 
21)EPIC-pndata. 

The EPIC-pn data are affected neither by telemetry drops 
nor by pile-up, in this observation, and can therefore be used 
for both spectral and timing analyses. However, as far as tim- 
ing/variability studies are concerned, some caution is necessary 
since burst-mode data suffer from the intrinsic photon loss re- 
lated to the 3% duty cycle (cp. Footnote 6). While the light 
curves only show tiny variations (up to a few percent, but com- 
patible with zero: see Fig. 2) at all timescales bigger than a 
few seconds, a quasi periodic oscillation (QPO) at about 0.6 
Hz, with a possibly harmonic signal at 1.2 Hz, is seen in the 
EPIC-pn data collected in burst mode (Fig. 3). This frequency 
is expected in x variability modes of GRS 1915+105, when 
the phenomenological correlation with the source flux is taken 
into account (cp. Fender & Belloni 2004). A forrest of narrow 
peaks at v > 100 Hz is apparent, due to spurious signals of 
instrumental origin. 

2.2.1. The RGS data 

Thanks to their excellent energy resolution, RGS cameras 
should be able to reveal subtle details in the 0.3-2.0 keV band. 
However, the source is so strongly absorbed below 1+2 keV 
that these data do not actually add much information. We used 
them mainly to check the EPIC-pn spectrum: in Fig. 4 we show 
the 0.5-2 keV spectra of both instruments, fitted with a sim- 
ple absorbed power law model, with F fixed to 2 (see Sect. 
3. 2). 8 The RGS and EPIC-pn spectra are clearly different, es- 
pecially with regard to the broad excess around 1 keV in the 
latter spectrum. The RGS spectrum shows a fast decline (more 
than two orders of magnitude, with very few counts at the lower 

8 Multiplicative constants to the RGS have also been added to allow 
for intercalibration problems (the best fit values for the constants are 
1.28 and 1.22 for the RGS1 and RGS2, respectively). RGS data, col- 
lected in the standard CCD read-out scheme, are marginally affected 
by pile-up, which amounts to a few percent and mainly affects the 
energies where the spectrum peaks. 




1 

Energy (keV) 



Fig. 4. The RGS and pn (burst mode) spectra fitted in the 0.5-2 keV 
range with a simple absorbed power law. The RGS and pn spectra 
clearly disagree each other. 



energies), without any apparent features. In Sect. 3.1.5 we go 
through the alternative hypotheses that can be invoked to ex- 
plain the RGS-pn discrepancy. 

RGS normally have a 15 ms timing resolution in HTR 
mode; we produced standard light curves, but did not get any 
useful results due to the low statistics on most of the band. 



3. Spectral analysis 

At the flux level of GRS 1915+105 (~ 10~ 8 erg cirT 2 s" 1 ), 
it is very likely that systematic errors cannot always be ne- 
glected with respect to the statistical ones. However, introduc- 
ing an energy-independent systematic uncertainty would pre- 
vent a meaningful use of the^ 2 statistics. Features at the 1-2% 
level must therefore be taken with caution as possibly of instru- 
mental origin (cp. Kirsch 2005). 

Even though the background issue is not relevant on most 
of the energy band, due to the extremely high S/N, nevertheless 
it may play a role at the highest energies (E > 10 keV), as well 
as at the energies where the source is strongly absorbed (E < 
1.5 keV). We therefore performed background subtraction; all 
reported results refer to background-subtracted analysis, unless 
stated otherwise. Of course, if a dust halo is present (cannot be 
checked with our data, given the lack of imaging capabilities 
of the timing modes; but see Greiner et al., 1998, on the issue), 
background subtraction may cause a spectral distortion in the 
soft band, depending on the distance from the source at which 
the background spectrum is taken (see below, Sect. 3.1.5). 

Spectra were rebinned in order to have about three bins 
per energy resolution element and at least 20 counts per bin at 
the same time. Spectral fits were performed with the XSPEC 
v.l 1.3.0 software package. In the following, all errors re- 
fer to a 90% confidence level for one interesting parameter 
(A;t 2 =2.71). 
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Fig. 5. The OB SI spectrum and data/model ratio when fitted with a 
simple power law and cold absorption with variable abundances. 




1 2 5 

Energy (keV) 



Fig. 6. The OBS 1 spectrum and data/model ratio for the best fit model 
(see text and Table 1). 



3.1. The timing-mode observation 

As discussed in the previous section, we are confident that the 
telemetry problems do not affect the spectrum significantly. 
Given the much larger total number of counts in OBS1, we 
decided to use it first for the spectral analysis. 

To illustrate the spectral complexity, in Fig. 5 we show the 
best fit spectrum and data/model ratio after fitting with a simple 
power law absorbed by cold matter. Element abundances were 
left free to vary with respect to the solar value ratios. Several 
features are apparent: an excess around 1 keV, several wiggles 
between ~ 1 .5 and 3 keV, an excess at the energies of the iron 
line complex, and a deficit above 8 keV. 

After several attempts, we found a - both statistically 
(^ 2 /d.o.f.=3 17.5/227) and physically - reasonable description 
of the spectrum in terms of (see Table 1 , and Fig. 6) a power law 
with r ~1.7, several absorption and emission features, an ion- 



ized and optically thick reflection component (and/or a warm 
absorber), and an optically thin reflection component to model 
the 1 keV excess. Given the fact that all remaining deviations 
of the data from the model are no larger than about 1-2%, i.e. 
on the order of the systematic errors in the calibration (which, 
for most of the energy range under analysis, is larger than the 
statistical errors), we did not try to improve the fit further. 

We now discuss the various spectral components separately 
for the sake of clarity. 

3.1 .1 . The continuum 

First of all, we note that a simple power law is basically suffi- 
cient to describe the continuum, apart from an additional com- 
ponent required to fit the >8 keV deficit and an apparent broad 
excess around ~ 1 keV. Substituting the power law with simple 
Comptonization models (like comptt or compst) does not lead 
to any significant improvement in the fit. No thermal emission 
from an accretion disk is apparent (the inclusion of a multicolor 
disk component is not required by the data). Of course, this 
strongly suggests that the innermost disk regions were absent 
during the observation. The apparent 1 keV excess is far too 
narrow to be fitted by even a single temperature blackbody (i.e. 
the narrower among physically motivated continuum models) 
or to be residual from "distorted" disk thermal emission. The 
nature of this feature will be discussed separately in Sect. 3. 1 .5. 

3.1.2. The cold absorber 

The photoelectric absorption has been fitted with the XSPEC 
model varabs, which assumes neutral matter and allows the 
elemental abundancies to vary separately. In practice, for the 
sake of simplicity, we grouped the elements so that within each 
group the elements are linked as: H, He, C, N, O; Ne, Na; Mg, 
Al; Si; S; CI, Ar, Ca, Cr; Fe, Co, Ni. The rationale behind this 
choice is both physical and practical. In fact, we grouped el- 
ements together whose origin is very likely common, but also 
those which are less abundant (and therefore cannot be mea- 
sured easily independently) with very abundant ones: e.g. Co 
and Ni with Fe. We also grouped H and He with C, N, and O 
because their effect is small in our energy range. 

In Table 1 the best fit values are reported. They are given in 
units of hydrogen equivalent column densities for solar abun- 
dances as in Anders & Grevesse (1989). A significant over- 
abundance of the heavier elements with respect to the lighter 
ones is apparent. These results are only partly consistent with 
those derived from the Chandra HETG observation (Lee et al. 
2002). In particular, Fe and S abundances are roughly the same, 
while our values for Mg are about twice, and 2/3 for Si. Fixing 
the abundances to the values of Lee et al. (2002), the fit is un- 
acceptable 0f 2 /d.o.f.=659. 9/231). We note that our observation 
has a lower energy resolution, but many more photons. 

These overabundancies with respect to the solar values 
(the differences are even greater assuming the ISM values, see 
Wilms et al. 2000) suggest that a significant fraction of the ab- 
sorber is local to the source. The value for the low Z elements 
(i.e. 1.6xl0 22 cirT 2 ) can then be assumed as the upper limit to 
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the ISM absorber column density. It is interesting to note that 
this value is similar to the one estimated by fitting the ~1 keV 
excess (see below). 

3.1.3. Reflection or warm absorption? 

The deficit above 8 keV can be fitted equally well by a simple 
smeared edge model (smedge) or by an ionized and relativis- 
tically blurred reflection component (refsch). As the smearing 
edge is not physically meaningful perse, but is usually assumed 
to represent the reflection component, we will not discuss the 
fits with this component further. 

The model adopted for the disk reflection component in- 
cludes the relativistic effects in Schwarzschild metric. The in- 
ner and outer radii are linked to those of the H-like iron line 
(see below). The disk inclination angle has been fixed to 70°, 
and the matter temperature to 10 6 K. The very large value of 
the inner radius (about 320r g ) justifies a posteriori the use of 
the Schwarzschild metric instead of the Kerr one (differences 
between the two metrics, even for a maximally rotating BH, 
are insignificant above ~10r g ), and it agrees with the lack of 
any detectable thermal disk emission. The iron abundance is 
five times the solar value, roughly consistent with the over- 
abundance found in the cold absorber (see above). The ion- 
ization parameter is pretty large, in agreement with the pres- 
ence of He- and H-like iron lines. The significant solid angle 
subtended by the disk to the illuminating source, i.e. 0.3x2^, 
implies that the emission occurs at distances comparable to the 
inner disk radius. The adopted model is necessarily a simplifi- 
cation of the real situation, as it adapts a unique ionization for 
the whole disk, while the ionization very likely decreases with 
increasing radius. The derived parameters, therefore, must be 
taken with caution and considered only as indicative. 

The presence of a H-like iron absorption line (see be- 
low) suggests warm absorption. Substituting the reflection 
component with a warm absorber (model absori), a signif- 
icantly worse, but not completely unacceptable fit is found 
(X 2 /d.o.f. =384.5/226), with a comparable iron overabundance 
(AFe=6.6) and power law photon index (T ~ 1.69). However, 
the column density (~ 7.7 x 10 21 cirT 2 ) is larger than that 
implied by the value of the EW of the absorption line (see 
Table 1), given the large iron overabundance (e.g. Bianchi et al. 
2005). Moreover, the presence of ionized disk lines (see below) 
requires the presence of the ionized reflection component. We 
thus fitted the spectrum with both the reflection and the warm 
absorber. No significant improvement is found with respect to 
the fit with the reflection only, but the upper limit to the column 
density, ~2xl0 21 cirT 2 , is now consistent with the absorption 
line EW. 

3.1 .4. The iron lines 

In Fig. 5, a clear broad excess at the energy of the iron line 
complex is apparent. The excess is at larger energies than 
6.4 keV, strongly suggesting ionized iron. We fitted this ex- 
cess with two unresolved emission lines {cr- 10 eV) peaking 
at 6.7 keV and 6.96 keV, which correspond to He- and H- 



like iron, respectively; and after Lee et al. (2002) we also in- 
cluded an unresolved absorption line with energy free to vary 
around the value of the H-like iron. The fit is unacceptable 
(y 2 /d.o.f.=617. 2/225), and not surprisingly, the absorption line 
was not found. A much better, but still hardly acceptable, fit 
(^ 2 /d.o.f.=425. 6/223) is found when leaving the widths of the 
lines free to vary. We found cr- 129 eV and 174 eV for the He- 
and H-like lines, respectively; an absorption line centered at 
~ 6.95 keV is also found. The high values of cr, along with 
the likely presence of the disk reflection component, strongly 
suggest that the lines are relativistically broadened: we there- 
fore substituted the broad Gaussians with diskline models, with 
the outer radius linked to that of the reflection component, as 
well as the inner radius of the H-like line (the inner radius of 
the He-like left free to vary, in the assumption that the matter is 
more ionized in the innermost regions). The results are reported 
in Table 1 . The fit is significantly better than with the Gaussians 
and - as expected on both physical grounds and from the fact 
that the width of the He-like line is smaller - the inner radius 
for this line is larger than that of the H-like line (the values 
of the radius must of course be taken merely as indicative, see 
discussion above). The energy of the narrow absorption line, 
6.95^J Q3 keV, is just below the rest frame value. 9 The line was 
assumed to be narrow (cr = 1 eV); when leaving cr free to vary, 
no improvement in the fit is found. 

We also tried to fit the iron emission with just one disk 
line model, with the energy free to vary. The fit is of compara- 
ble statistical quality, but the line energy, 6.86+0.01 keV, does 
not correspond to any rest frame atomic energy of significance. 
Moreover, the good fit is obtained at the expense of a very large 
value for the iron abundance of the reflection component (~18 
times the solar value). Therefore, the model with two disk lines 
has to be preferred. 

The EW of both lines (28 eV) seems rather small when 
compared to the amount of the reflection component (Matt et al. 
1996). However, one must note again that the ionization most 
likely changes along the radius in the disk. While reflection 
occurs (even if with changing spectra) throughout the disk, line 
emission occurs only if the corresponding ion is present, which 
probably only happens in a fraction of the disk. 

3.1 .5. The ~1 keV excess 

In Sect. 2.2.1 we have shown that, during the burst-mode ob- 
servation, the RGS and pn disagree, with the former not show- 
ing the 1 keV excess. (Unfortunately, during the timing-mode 
observation the RGS were switched off). There are several pos- 
sible explanations for both the excess and the RGS-pn discrep- 
ancy. As we will see, none of them is fully satisfactory, and we 
consider the issue still open as we wait for improvements in 
calibrations. 

A pn instrumental origin for the broad excess seems un- 
likely for various reasons. The amplitude (about 40% of the un- 
derlying continuum, see the data vs. model ratio in Fig. 5) is far 
too large, given the known calibration uncertainties (cp. Kirsch 

9 Actually, this is a blend of two lines, with centroid energy at 6.966 
keV. 
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2005). Moreover, the excess is present in the Timing as well as 
in the burst mode, and was also observed with BeppoSAX (cp. 
Fig. 3 in: Feroci et al. 1998). A similar feature has also been 
seen in other binaries (see e.g. the case of the dipping LMXB 
4U 1323-62: Boirin et al. 2005), but not in all; in particular, it 
is not present in the burst mode observation of the Crab Nebula 
(Kirsch et al. 2005). 

However, given the unusual brightness of our source and 
the fact that the excess arises in the part of the spectrum 
where the count rate drops precipitously due to absorption, a 
calibration problem that is possibly related to the redistribution 
matrix cannot be completely excluded at the moment. 

Let us then assume that this excess is real. We could not fit 
this feature with any continuum model put behind the whole 
absorber. As seen in Sect. 3.1.1, it cannot be fitted with black- 
body models; the only model able to fit it is a broad Gaussian 
with an EW of several tens of keV, which has no physical expla- 
nation. Good fits are instead found if the emission is screened 
by a thinner absorber. There are three possibilities in this re- 
spect: (i) a confusing source, (ii) a dust scattering halo, or (iii) 
a component related to the system but lying outside the intrin- 
sic absorber. 

The confusing source hypothesis seems unlikely. In tim- 
ing modes, source data are extracted from a narrow strip (here 
about 16"); the required, "polluting" source would have to be 
bright (0.5-2 keV observed flux of about 4xl0~ 12 erg crrT 2 
s _1 ): but no source is present in the ROSAT All Sky Survey 
within a radius of 0.3° around GRS 1915+105 (we checked 
with the HEASARC tool PIMMS that such a source would be 
well above the typical RASS detection limit; of course, the pos- 
sibility of a variable source cannot be ruled out). 

Since the dust scattering is energy dependent, ignoring a 
possible halo component leads to a distortion of the spectrum in 
a small energy band, because of an incorrect background sub- 
traction (depending on the distance from the source at which 
the background spectrum is taken); however, less soft photons 
would be "observed". A possibility is that we are including the 
halo in the extraction region. Since there is no obvious dust 
scattering model to fit the data, imaging may be the only way 
to test the halo hypothesis: given the particular modes adopted, 
this sort of analysis is very difficult (and altogether impossible 
in the burst mode). We only note that the dust scattering halo 
observed in Cyg X-2 (Costantini et al. 2005) has quite a broad 
spectrum (when normalized to the source spectrum), broader 
than the feature under discussion. Another problem with the 
scattering halo hypothesis is that the relative contribution of 
the excess is similar in the timing- and burst-mode observa- 
tions, which is possible only if most of the halo is within 16", 
i.e. the width of the strip from which the timing-mode spectrum 
was extracted. 

Even if unsatisfactory in many respects, the two above- 
mentioned hypotheses have the merit of invoking extended 
or off-axis emission, which can in principle help explain the 
presence of the excess in the pn spectra but not in the RGS. 
However, as the RGS spectrum is extracted over a region of 
about 2.5 arcminutes, the halo hypothesis seems to be excluded 
for the reasons discussed above. A confusing source, provided 



it lies at a greater distance, is instead still tenable. We note 
that all Chandra observations of GRS 1915+105 have been 
performed in the CC (i.e. non imaging) mode, therefore the 
above-mentioned hypotheses could not be tested yet. 

In the assumption that the RGS spectrum suffers from in- 
strumental/calibration problems, the 1 keV component could 
be related to the system, and a possible candidate would be the 
emission by a disk wind. Both a collisionally ionized plasma 
model and a photoionized plasma model do indeed fit the 
EPIC-pn data well. 

In the first case, the temperature of the plasma is about 1.3 
keV (obtained using the mekal model), with metal abundance 
1 .2 times the solar value; the ISM absorber has an /V H of 9x 10 21 
crrT 2 . The emission integral is 7.2xl0 57 <f 2 crrT 3 (where d\o is 
the distance to the source in units of 10 kpc), which, assum- 
ing that the density of the matter is constant along the emitting 
region, translates into a mass of the wind M w x 6 M & d\ { jT x , 
where n is the number density. Assuming that the size of the 
disk wind is at most on the order of the binary separation (about 
5xl0 10 cm according the binary system parameters given by 
Harlaftis & Greiner 2004), we obtain a lower limit to n of about 
8xl0 12 <fio crrT 3 . With this size, the Thomson optical depth 
would already be substantial (~ 0.25), and be increasing fur- 
ther with decreasing size. Thus, the matter would not be really 
optically thin, so this solution is not self-consistent. 

Alternatively, the matter may be in photoionization equilib- 
rium with the emission due to reflection of the primary X-ray 
radiation. To model this case, we adopted a power law with 
the same photon index as the primary one (to model Thomson 
scattered radiation) plus a Gaussian line. In this case, the ISM 
Nu is larger, ~1.5xl0 22 crrT 2 (this value may be assumed as 
an estimate of the ISM column density) 10 . From the observed 
flux, the optical depth of the wind is about 0.001//, where / 
is the covering factor of the wind to the primary radiation. The 
centroid energy of the Gaussian line (~0.97 keV), its width (90 
eV), and its EW against the reflected continuum (5.6 keV), all 
suggest a blend of Ne K and Fe L lines. 

3.1.6. Other emission/absorption lines 

Other emission and absoption lines are required to fit the spec- 
trum. They are: two unresolved emission lines at about 1.85 
keV and 2.25 keV, with EW=19 eV and 8 eV, respectively; two 
broad absorption lines at about 2.5 and 2.93 ke V, with <x of 0. 1 1 
and 0.14 keV and EWs of -17 and -30 eV, respectively. 

Some of these features may be instrumental artifacts, since 
this energy range is notoriously troublesome due to the changes 
in the effective area. 11 However, some of them appear too 

10 For the sake of simplicity we fitted the absorption with a pure 
gas photoelectric absorption (model wabs in XSPEC). Using a more 
appropriate ISM model that includes the effects of dust (TBabs in 
XSPEC, Wilms et al. 2000) does not change the results significantly. 
The fit returns basically the same x 1 and the same column density. 

1 1 See Figs.30 and 3 1 in the XMM-Newton User's Handbook: http:// 
xmm.vilspa.esa.es/external/xmm_user_support/documentation/ 
uhb/node32.html . 
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prominent to be entirely due to calibration problems. The first 
emission line can be readily identified with the Si xm recom- 
bination line. If arising from the same photoionized material 
responsible for the ~1 keV excess (see above), its total EW 
would be about 16 keV, a very large but not fully unrealistic 
value, given the large overabundance we found in the local ab- 
sorber and provided that the reflecting matter is still optically 
thin to resonant scattering (Matt et al. 1996), which is possible 
if the matter is largely turbulent. The same is true for the sec- 
ond emission line, which can be identified with Si xiv and S xn, 
even if the best fit value of the energy is slightly lower. 

More difficult to explain are the two broad absorption lines. 
The energies correspond to the Ka and K/3 recombination lines 
from S xv. In this case, the absorbing material cannot be the 
same as that responsible for the Fe xxvi absorption line (see 
Sect. 3.1.4), since the latter is more ionized. It is also worth 
recalling that the modeling of the absorption edges is oversim- 
plified, because it does not take into account X-ray absorption 
fine structures (XAFS), which are generated if the material is 
in grains (e.g. Lee et al. 2002); this may introduce spurious 
features in the residuals, even if it cannot fully explain them. 
Moreover, we again stress that it is not inconceivable that a 
significant part of these features are of instrumental origin. 
Data with higher energy resolution and optimized calibration 
are necessary to solve this issue. 

3.2. The burst-mode observation 

We fitted the burst-mode spectrum with the same model used 
for the timing-mode spectrum, fixing those parameters that are 
not expected to vary. In particular, we fixed the interstellar ab- 
sorption column density, the ratio between the elemental abun- 
dancies of the local absorber, and the iron abundance of the 
reflection component at the best fit values obtained for the 
timing-mode observation. A good fit (y 2 /d.o.f.=248. 2/219) is 
found. The best fit parameters are summarized in Table 2. 

In comparison with the timing-mode observation (see 
Fig. 7), the spectrum is steeper (T ~2 instead of ~1.7), and 
the reflection component much stronger (R ~ 1.7 instead of 
0.3) and more ionized (£ ~ 3300 instead of 940 erg cm s~'). 
The inner and outer radii are much smaller (r, <20, r Q ~100). 
The local absorber is slightly larger, while the ionized iron disk 
lines are not detected (upper limits to the EW of about 10 eV for 
both the He- and H-like lines. The Fe H-like absorption line 
is still present, with a centroid energy of 6.98±0.02, i.e. consis- 
tent with both the rest frame value and the value obtained in the 
timing-mode observation. The EW is larger, about -50 eV, due 
to the fact that the line is now best fitted with a resolved width 
cr of 0. 1 1 +0.03 keV. The 1 .85 and 2.24 keV emission lines are 
not required, with upper limits to the EWs of 5 and 2 eV, re- 
spectively. A 6.4 keV iron emission line is instead marginally 
found, with an EW of 6+4 eV. 

The above-mentioned results are consistent with a picture 
in which the disk is now more extended downwards and ionized 
(in fact, so ionized that no significant relativistic iron line is ob- 
served). The steeper spectrum can be due to the more efficient 
cooling, itself due to the increase in soft disk photons. The fact 
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Fig. 7. The spectrum and data/model ratio of the burst-mode observa- 
tion with the best fit model for the timing-mode observation. Clearly, 
iron line emission is over-subtracted, and the deficit at higher energies 
turns out to be more important than in OBS1; note again the strong 
excess at ~ 1 keV and the wiggles at ~ 2-3 keV (cp. with Fig. 5). 

that no thermal disk emission is required, however, implies that 
the values for the inner and outer radii must be used with cau- 
tion. It must be noted that they are now determined only via the 
reflection component, with all the associated uncertainties (see 
discussion in the previous section), as no disk line is present. 

However, these results, especially those regarding the re- 
flection components, can at least partly be an artifact due to 
calibration problems. A deficit of photons above 6+7 keV, i.e. 
similar to what we find (see Fig. 7), is indeed also present in 
the burst-mode observation of the Crab Nebula (Kirsch et al. 
2005), a source where of course no significant reflection com- 
ponent is expected. 

The 1 ke V excess, which still persists, can again be well fit- 
ted by a power law plus an emission line, the latter with a flux 
slightly larger than that obtained in the timing-mode observa- 
tion. In the reflection scenario discussed in the previous sec- 
tion, the optical depth is about a factor 10 larger, which implies 
a much lower EW of the line. These findings can be explained, 
at least qualitatively, by assuming that the matter is now thicker 
to the lines, with consequent self-absorption effects (see e.g. 
Matt et al. 1996). 

4. Conclusions 

In this paper for the first time we have presented XMM-Newton 
observations of GRS 1915+105. Several problems make it dif- 
ficult to observe the source fruitfully with this satellite, first of 
all due to the extremely high flux - in our case: F2-10 keV ~ 
6.0 + 6.6 (unabsorbed: ~ 8.7 + 10.7) xl0~ 9 erg cirT 2 s" 1 - and 
to the notorious, erratic source variability, which does not help 
in triggering observations. We succedeed at both a) observing 
the source in a well-defined, stable physical/spectral state and 
b) collecting EPIC-pn useful data, only marginally corrupted 
by telemetry problems. 
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In both observations (April 17th and 21st, 2004), the source 
has been caught in a "plateau" state, also characterized by 
stronger jet activity and higher radio level, which we identify 
with the conventional "C" state j x variability mode as defined 
by Belloni et al. (1997a,b, 2000). In this state the source shows 
a QPO at ~ 0.6 Hz, i.e. what is expected in plateau intervals 
when the frequency vs. spectral hardness correlation is taken 
into account. 

In order to try to disentagle the different spectral compo- 
nents, and thus to better discriminate among the viable spec- 
tral/physical models, we used the rms vs. E method by Ponti 
et al. (2004). The resulting rms is lower than 0.1 across the en- 
ergy band. We further analyzed the variable F var , which takes 
the Poisson noise into account, and concluded that this variable 
is compatible with being null on the whole energy band, i.e. all 
spectral components are compatible with being constant during 
the two observations. 

A priori, some of the features in the spectrum may be af- 
fected by dust halo scattering, although this cannot explain the 
1 keV excess (see Sect. 3.1.5). However, testing this hypothesis 
would first of all require an accurate analysis based on proper 
source imaging (if available); further, spectral modelling of 
such effects is not easy (see e.g. Costantini et al. 2005). Being 
beyond the scope of this paper, this study should be addressed 
in future works. 

We adopted a power law continuum model. This could be 
related to emission by a hot corona or to Comptonized thermal 
emission e.g. from the jet basis (as proposed by Rodriguez et 
al., 2004); however, an optically thick reflector is required to 
account for the smeared edge at ~ 7 keV (with a covering ra- 
tio of ~ 0.4 -r 1.7 x 2n in the two observations, respectively). 
This may be evidence of an accretion disk being present, or just 
optically thick, only at quite a large distance from the central 
compact object, at least in the first observation (r;/r g > 300 
in OBS1, ~ 20 in OBS2). The relatively large amount of the 
reflection components implies that the primary X-ray emitting 
region should have a size that is at least comparable to the in- 
ner disk radius. That the disk is truncated, i.e. not present in 
the innermost part, is also suggested by the non-detection of a 
thermal disk emission component. 

Several line residuals are superimposed on the modeled 
continuum. Part of these may be due to calibration uncertain- 
ties, especially at the energies where changes in the EPIC ef- 
fective area take place (e.g. 1-3 keV). However, we found clear 
evidence of ionized iron emission around ~ 7 keV: data are well 
fitted with two ionized Fe Ka lines, possibly affected by mild 
relativistic broadening (being produced far away from the BH 
event horizon), plus a narrow absorption feature at ~ 6.95 keV. 

The puzzling presence of an intense, broad excess around 1 
keV can be explained in terms of reflection by an optically thin 
wind, provided that the feature is real and that the discrepancies 
between the RGS and pn spectra (discussed in Sections 2.2. 1 
and 3.1.5) have to be attributed to RGS calibration problems. 
This assumption is still unsubstantied, and should be regarded 
as a mere working hypothesis. The excess is indeed well fitted 
with a power law plus a line, unobscured by material that is 
intrinsic to the system. The centroid energy of the Gaussian 
line (-0.97 keV), its width (90 eV), and its EW against the 



reflected continuum (5.6 keV), together point to a blend of Ne 
K and Fe L lines. The value of the equivalent H column density 
is interestingly similar to the value of the obscuration by low Z 
elements (H, He, C, N, O) at the source core -N H ~ 1 .6 X 10 22 
cm -2 : in the disk wind hypothesis, this may therefore be taken 
as an upper limit to the interstellar matter column density. This 
value matches the expected galactic absorption in that direction 
well (Dickey & Lockman 1990). 

On the other hand, a significant fraction of the absorber 
must be local to the source. We adopted a variable absorption 
model (varabs in XSPEC), assuming neutral matter and group- 
ing the elements on the basis of both physical and practical con- 
siderations: elements which have probably a common origin, 
but also elements which are not very abundant (and therefore 
cannot be easily measured independently one from the other) 
with very abundant ones (e.g. Co and Ni with Fe): see the re- 
sults in Table 1 . A significant overabundance of the heavier el- 
ements compared to the lighter ones is apparent, which sug- 
gests that a significant fraction of the absorber, traced by heav- 
ier species, is local to the source. These results are somewhat 
different from those derived from Chandra HETG data (Lee et 
al. 2002): in particular, Fe and S abundances are roughly the 
same, while our values for Mg are about twice, and 2/3 for Si. 
The overall scenario is, however, very similar. 

Clearly, the intrinsic absorption may be subject to substan- 
tial changes on longer timescales, as already observed with 
Rossi-XTE in similar plateaux (class x, state C): from a mean 
value of ~ 4.5 x 10 22 cirT 2 , the total A^h may occasionally rise 
up to ~ 6 x 10 22 crrT 2 (Belloni et al. 2000). 

Assuming the customary value of ~ 12.5 kpc for the 
source's distance (Rodriguez & Mirabel, 1999) we obtain an 
estimate of the intrinsic luminosity ^"fokev = 1 .7 -r 2.0 X 10 38 
erg s _1 in the two observations. 12 
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N, 



H,He,C,N,0 
A'Mg.Al 

N si 
N s 



N, 



Cl,Ar,Ca,Cr 



Fe.Co.Ni 



Ei 
F, 
EW 

E, 
F, 
EW 



Ei 
F, 
EW 
cr 



n/r g 

F, 
EW 



E c 

n/r 8 

F, 
EW 



R/2n 

n/r g 

F2-IO keV 

X 2 /d.o.f. 



Power law 

1.686^ 

Cold absorption 

I. 60 + ^x 10 22 cm- 2 
7.46 + ^i x 10 22 cm- 2 
7.57 + ^ x 10 22 cm- 2 

5 - 70 -oi2 x 1022 cm " 2 
4.69$% x 10 22 

II. + !'ixl0 22 



-2 



1 1.7+°f xl0 2; 



cm 



Emis. lines 

L846 -aoM keV 

2 ph cm -2 : 



1.35^ x 10 



19 eV 



2.244!°,;™ keV 



0.010 

3.62:»;« x 10- 3 
8eV 



ph cm 



Abs. line 

6.95 + ^ keV 
-7.9+0- 2 x 10 J4 ph cm- 2 s" 1 
-9eV 
1 eV (frozen) 

1st disk line (He-like Fe) 

6.7 keV (fixed) 

2.37:°-;* x 10- 3 ph cm- 2 s" 1 
28 eV 

2nd disk line (H-like Fe) 

6.96 keV (fixed) 

320Z 
2.20$% x 10- 3 



ph cm 
28 eV 



Reflection 

35 +a02 

5 2 +av 

940!™ erg cm s _1 
320+so 

-^"-60 



0.6 (unabs: 



0.87) 10- 8 erg cm" 2 s" 1 



317.5/227 



r 


Power law 

2 04 +0 01 

-0.02 


^H,He,C,N,0 


Cold absorption 

1.98^02 x 10 cm 


R/2n 

r i/r g 


Reflection 

1 69 +016 
<20 
93 + |^ 
3300!™° erg cm s" 1 


F2-W keV 


~ 6.6 (unabs: ~ 10.7) 10 -9 erg cm" 2 s" 1 


^/d.O.f. 


248.2/219 



Table 2. The best fit model parameters for the burst-mode observation 
(OBS2). See the text (Sect. 3.2) for details. 



Table 1. The best fit model parameters for the timing-mode observa- 
tions. The column densities are the Hydrogen equivalent column for 
solar abundances as in Anders & Grevesse (1989). The outer radii of 
the reflection and disk line models have been linked together in the fit, 
yielding r a /r g = 900!™; r, in the 2nd disk line model has been linked 
together to r, of the reflection component. See the text for further dis- 
cussion. 



